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We compare waveforms obtained by numerically evolving nonspinning binary black holes to post- 
Newtonian (PN) template families currently used in the search for gravitational waves by ground- 
based detectors. We find that the time-domain 3.5PN template family, which includes the inspiral 
phase, has fitting factors (FFs) > 0.96 for binary systems with total mass M = 1O-2OM0. The time- 
domain 3.5PN effective-one-body template family, which includes the inspiral, merger and ring-down 
phases, gives satisfactory signal-matching performance with FFs > 0.96 for binary systems with total 
mass M = 1O-12OM0. If we introduce a cutoff frequency properly adjusted to the final black-hole 
ring-down frequency, we find that the frequency-domain stationary-phase-approximated template 
family at 3.5PN order has FFs > 0.96 for binary systems with total mass M = 10-20Mq. However, 
to obtain high matching performances for larger binary masses, we need to either extend this family 
to unphysical regions of the parameter space or introduce a 4PN order coefficient in the frequency- 
domain GW phase. Finally, we find that the phenomenological Buonanno-Chen-Vallisneri family 
has FFs > 0.97 with total mass Af = 10-120Mq. The main analyses use the noise spectral-density 
of LIGO, but several tests are extended to VIRGO and advanced LIGO noise-spectral densities. 

PACS numbers: 04.25.Dm, 04.30.Db, 04.70.Bw, x04.25.Nx, 04.30.-w, 04.80.Nn, 95.55.Ym 
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I. INTRODUCTION 

The search for gravitational-waves (GWs) from coa- 
lescing binary systems with laser interferometer GW de- 
tectors [l|, [1, H, 0, § is based on the matched-filtering 
technique, which requires accurate knowledge of the 
waveform of the incoming signal. In the last couple of 
years there have been several breakthroughs in numeri- 
cal relativity (NR) [1,0, S], and now independent groups 
are able to simulate the inspiral, merger and ring-down 
phases of generic black- hole (BH) merger scenarios, in- 
cluding different spin orientations and mass ratios 
However, the high computational cost of running such 
simulations makes it difficult to generate sufficiently long 
inspiral waveforms that cover the parameter space of as- 
trophysical interest. 

References [HI found good agreement between ana- 
lytic (based on the post-Newtonian (PN) expansion) and 
numerical waveforms emitted during the inspiral, merger 
and ring-down phases of equal-mass, nonspinning binary 
BHs. Notably, the best agreement is obtained with 3PN 
or 3.5PN adiabatic waveforms [l^l (henceforth denoted 
as Taylor PN waveforrns) and 3.5PN effective-one-body 
(EOB) waveforms [H, M, M, M, [13, [3, [ll- In addi- 
tion to the inspiral phase the latter waveforms include 
the merger and ring-down phases. Those comparisons 
suggested that it should be possible to design hybrid nu- 
merical/analytic templates, or even purely analytic tem- 
plates with the full numerics used to guide the patching 
together of the inspiral and ring-down waveforms. This is 



an important avenue to template construction as eventu- 
ally thousands of waveform templates may be needed to 
extract the signal from the noise, an impossible demand 
for NR alone. Once available, those templates could be 
used by ground-based laser interferometer GW detectors, 
such as LIGO, VIRGO, GEO and TAMA, and in the fu- 
ture by the laser interferometer space antenna (LISA) for 
detecting GWs emitted by solar mass and supermassive 
binary BHs, respectively. 

This paper presents a first attempt at investigating the 
closeness of the template families currently used in GW 
inspiral searches to waveforms generated by NR simula- 
tions. Based on this investigation, we shall propose ad- 
justments to the templates so that they include merger 
and ring-down phases. In contrast, Ref. 21i] examined 
the use of numerical waveforms in inspiral searches, and 
compared numerical waveforms to the ring-down tem- 
plates currently used in burst searches. Similar to the 
methodology presented here, fitting factors (FFs) [see 
Eq. ([2]) below] are used in Ref. [2l| to quantify the ac- 
curacy of numerical waveforms for the purpose of de- 
tection, as well as the overlap of burst templates with 
the waveforms. Reference [2l| found that by computing 
FFs between numerical waveforms from different reso- 
lution simulations of a given event, one can recast the 
numerical error as a maximum FF that the numerical 
waveform can resolve. In other words, any other tem- 
plate or putative signal convolved with the highest reso- 
lution numerical simulation that gives a FF equal to or 
larger than this maximum FF is, for the purpose of de- 
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tection, indistinguishable from the numerical waveform. 
We will explore this aspect of the problem briefly. The 
primary conclusions we will draw from the analysis do 
not crucially depend on the exactness of the numerical 
waveforms. What counts here is that the templates can 
capture the dominant spectral characteristics of the true 
waveform. 

For our analysis we shall focus on two nonspinning 
equal-mass binary simulation waveforms which differ in 
length, initial conditions, and the evolution codes used 
to compute them: Cook-Pfeiffer quasi-eq uilibrium ini- 
tial data built on Rcfs. [H, [H, [H, [H, evolved with 
Prctorius' generalized harmonic code 6|, and Brandt- 
Briigmann puncture data [13] evolved using the Goddard 
group's moving-puncture code Q. We also consider two 
nonspinning unequal-mass binary simulations with mass 
ratios 7712/"^! ~ 1-5 and m2/mi = 2 produced by the 
Goddard group. 

The paper is organized as follows. In Sec. |TT] we dis- 
cuss the phase differences between PN inspiraling tem- 
plates. In Sec. mil we build hybrid waveforms by stitch- 
ing together PN and NR waveforms. We try to under- 
stand how many NR cycles are needed to obtain good 
agreement between NR and PN waveforms, to offer a 
guide for how long PN waveforms can be used as accu- 
rate templates. In Sec. IIVI we compute the FFs between 
several PN template families and NR waveforms. We 
first focus on low-mass binary systems with total mass 
M = 1O-3OM0, then high-mass binary systems with to- 
tal mass M = 3O-I2OM0. Finally, Sec. |V] contains our 
main conclusions. In Appendix [X] we comment on how 
different representations of the energy-balance equations 
give GW frequencies closer to or farther from the NR 
ones. 



II. PHASE DIFFERENCES IN 
POST-NEWTONIAN INSPIRALING MODELS 



Starting from Ref. [28|, which pointed out the impor- 
tance of predicting GW phasing with the highest possi- 
ble accuracy when building GW t emp lates, many subse- 
quent studies [13, M, M, M, m, [M m, [3| (those ref- 
erences are restricted to the nonspinning case) focused 
on this issue and thoroughly tested the accuracy of those 
templates, proposing improved representations of them. 
These questions were motivated by the observation that 
comparable-mass binary systems with total mass higher 
than SOA^Q merge in-band with the highest signal-to- 
noise ratio (SNR) for LIGO detectors. It follows that the 
corresponding templates demand an improved analysis. 

In the absence of NR results and under the urgency 
of providing templates to search for comparable-mass bi- 
nary BHs, the analytic PN community pushed PN cal- 
culations to higher PN orders, notably 3.5PN order p^ . 
and also proposed ways of resumming the PN expan- 
sion, either for conservative dynamics (the FOB ap- 
proach [1^, [3, [13] ) , radiation-reaction effects (the Padc 



resummation [lEl), or both [Tj,[Tal- Those results lead to 
several conclusions: (i) 3PN terms improve the compar- 
ison between analytic and (numerical) quasi-equilibrium 
predictions [2^, [26|, [13, [IH] ; (ii) Taylor expanded and re- 
summed PN predictions for equal-mass binary systems 
are much closer at 3.5PN order than at previous PN 
orders, indicating a convergence between the different 
schemes [TSi . i2ft . i31i] ; (iii) the two-body motion is quasi- 
circular until the end of a rather blurred plunge [l3| , (iv) 
the transition to ring-down can be described by an ex- 
tremely short merger phase [l3 . [Tsj . Today, with the NR 
results we are in a position to sharpen the above con- 
clusions, and to start to assess the closeness of analytic 
templates to numerical waveforms. 

Henceforth, we restrict the analysis to the three time- 
domain physical template families which are closest to 
NR results [13, [Ol: the adiabatic Taylor PN model 
(Tpn) [see, e.g., Eqs. (1), (10), and (11)-(13) in Ref. ^] 
computed at 3PN and 3.5PN order, and the nonadia- 
batic FOB model (Epn) [see e.g., Fqs. (3.41)-(3.44) in 
Ref. [l3|] computed at 3.5 PN order. We shall denote 
our models as Tpn(n) and Fpn(n), n being the PN or- 
der. The Tpn model is obtained by solving a particular 
representation of the balance equation. In Appendix |^ 
we briefly discuss how time-domain PN models based on 
different representations of the energy-balance equation 
would compare with NR results. 

The waveforms we use are always derived in the so- 
called restricted approximation which uses the amplitude 
at Newtonian order and the phase at the highest PN 
order available. They are computed by solving PN dy- 
namical equations providing the instantaneous frequency 
Lo{t) and phase (p{t) ^ (po + j^^ uj{t')dt' , thus 



h{t) = ytcj(i)2/3 cos[20(t)] 



(1) 



where to and 0o are the initial time and phase, respec- 
tively, and ,4 is a constant amplitude, irrelevant to our 
discussion. The inclusion of higher-order PN corrections 
to the amplitude can be rather important for certain 
unequal-mass binary systems, and will be the subject of 
a future study. 

When measuring the differences between waveforms 
wc weight them by the power spectral-density (PSDs) of 
the detector, and compute the widely used fitting factor 
(FF) (i.e., the ambiguity function or normalized over- 
lap), or equivalently the mismatch defined as 1-FF. Fol- 
lowing the standard formalism of matched-filtering [see, 
e.g., Refs. [13, dll, [3^], we define the FF as the overlap 
(/ii(i), /i2(t)) between the waveforms hi{t) and h2{t): 

(/.i(t),/.2(t))^4Re TMP^^d/, 



FF = max 



Jo Snif) 
(/ii,/i2(io,0o, A')) 



*o.<^o,A* y/{huhi){h2{to, (Ao, A»), h2{to, (j)o, X')) 



(2) 



where hi{f) is the Fourier transform of hi{t), and Sh{f) 
is the detector's PSD. Thus, the FF is the normalized 
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overlap between a target waveform hi (t) and a set of tem- 
plate waveforms /i2(^0j 4*0 ^ maximized over the initial 
time to, initial phase t^o, and other parameters A*. Some- 
times we are interested in FFs that are optimized only 
over to and (j>o; we shall denote these as FFq. For data- 
analysis purposes, the FF has more direct meaning than 
the phase evolution of the waveforms, since it takes into 
account the PSDs and is proportional to the SNR of the 
filtered signal. Since the event rate is proportional to 
the cube of the SNR, and thus to the cube of the FF, a 
FF= 0.97 corresponds to a loss of event rates of ^ 10%. 
A template waveform is considered a satisfactory repre- 
sentation of the target waveform when the FF is larger 
than 0.97. 

When comparing two families of waveforms, the FF is 
optimized over the initial phase of the template wave- 
form, and we also need to specify the initial phase of the 
target waveform. Since there is no preferred initial phase 
of the target, two options are usually adopted: (i) the ini- 
tial phase maximizes the FF or (ii) it minimizes the FF. 
The resul ting FFs are referred to as the best and mini- 
max FFs [29| . All FFs we present in this paper are min- 
imax FFs. Although the FF of two waveform families is 
generally asymmetric under interchange of the template 
family [3l| . the best and the minimax FFqS are symmet- 
ric (see Appendix B of Ref . [1^ for details). Henceforth, 
when comparing two waveform families using FFq , we do 
not need to specify which family is the target. 

We shall consider three interferometric GW detectors: 
LIGO, advanced LIGO and VIRGO. The latter two have 
better low-frequency sensitivity and broader bandwidth. 
For LIGO, we use the analytic fit to the LIGO de- 
sign PSD given in Ref. [g^l; for advanced LIGO we use 
the broadband configuration PSD given in Ref. [13]; for 
VIRGO we use the PSD given in Ref. 

In Fig. [TJ we show the FFqS as functions of the ac- 
cumulated difference in the number of GW cycles be- 
tween waveforms generated with different inspiraling PN 
models and for binary systems with different component 
masses. We first generate two waveforms by evolving 
two PN models, say, "PNi" and "PN2" which start at 
the same GW frequency /gw — BOHz and have the same 
initial phase. The two waveforms are terminated at the 
same ending frequency /gw = /end up to a maximum 

/end, max = min(/ond,PNi , /end.PNs ) ; wllCrC /cnd,PN IS the 

frequency at which the PN inspiraling model ends. (For 
Tpn models this is the frequency at which the PN energy 
has a minimum; for Epn models it is the FOB light-ring 
frequency.) Then, we compute the difference in phase 
and number of GW cycles accumulated until the ending 
frequency 

ANgW = — = - [0PNi(/cnd) - 0PN2(/ond)] ■ (3) 
TT TT 

By varying /end (up to /end, max) ANgw changes, though 
not necessarily monotonically. Although there seems to 
be a loose correlation between the FFqS and ANgw; it 
is hard to quantify it as a one-to-one correspondence. 
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FIG. 1: We show FFqs between waveforms generated from 
the three PN models Tpn(3), Tpn(3.5) and Epn(3.5) versus 
ANgw [see Eq. Q]. The FFqS are evaluated with LIGO's 
PSD. Note that for Tpn(3.5) and Epn(3.5) and a {15 + 3)Mq 
binary, the lowest FFq is 0.78 and the difference in the number 
of GW cycles ANqw - 2. In the limit ANgw ^ 0, the FFq 
goes to unity. 



For example, a phase difference of about half a GW cy- 
cle (ANgw — 0.5) is usually thought to be significant. 
However, here we find relatively high FFqS between 0.97 
and > 0.99, depending on the masses of the binary and 
the specific PN model used. This happens because the 
FF between two waveforms is not determined by the to- 
tal phase difference accumulated, but rather by how the 
phase difference accumulates across the detector's most 
sensitive frequency band. The relation between FFs and 
phase differences is also blurred by the maximization 
over the initial time and phase: shifting the phase by 
half a cycle from the most sensitive band to a less sen- 
sitive band can increase the matching significantly. We 
conclude that with LIGO's PSD, after maximizing only 
on initial phase and time, Epn(3.5) and Tpn(3.5) tem- 
plates are close to each other for comparable-mass bi- 
nary systems M = 6-30Mq with FFq > 0.97, but they 
can be different for mass ratios m2/mi ~ 0.3 with FFq 
as low as ~ 0.8. Tpn(3) and Tpn(3.5) templates have 
FFq ^0.97 for the binary masses considered. Note that 
for m.2/mi = 1 [~ 0.3] binary systems, Tpn(3.5) is closer 
to Epn(3.5) [Tpn(3)] than to Tpn(3) [Epn(3.5)]. Note 
also that when maximizing on binary masses the FFs can 
increase significantly, for instance, for a (15 -I- 3)Mq bi- 
nary, the FF between Tpn(3.5) and Epn(3.5) waveforms 
becomes > 0.995, whereas FFq ~ 0.8. 



III. BUILDING AND COMPARING HYBRID 
WAVEFORMS 

Recent comparisons [lol [Tl| between analytic and nu- 
merical inspiraling waveforms of nonspinning, equal-mass 
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binary systems have shown that numerical waveforms are 
in good agreement with Epn(3.5), Tpn(3) and Tpn(3.5) 
waveforms. Those results were assessed using eight and 
sixteen numerical inspiral GW cycles. Can we conclude 
from these analyses that Epn(3.5), Tpn(3.5) and Tpn(3) 
can safely be used to build a template bank for detecting 
inspiraling GW signals? A way to address this question 
is to compute the mismatch between hybrid waveforms 
built by attaching either Epn or Tpn waveforms to the 
same numerical waveform, and varying the time when 
the attachment is made. This is equivalent to varying 
the number of numerical GW cycles n in the hybrid tem- 
plate. The larger n the smaller the mismatch, as we are 
using the same numerical segment in both waveforms. 
For a desired maximum mismatch, say 3%, we can then 
find the smallest number n of numerical cycles that is 
required in the hybrid waveform. This number will, of 
course, depend on the binary mass and the PSD of each 
detector. 



A. Hybrid waveforms 

We build hybrid waveforms by connecting PN wave- 
forms to NR waveforms at a chosen point in the late 
inspiral stage. As mentioned before, we use NR wave- 
forms generated with Pretorius' code and the God- 
dard group's [s^ code. Pretorius' waveform is from an 
equal-mass binary with total mass M, and equal, co- 
rotating spins (a = 0.06). The simulation lasts ~ 671M, 
and the waveform has ~ 8 cycles before the formation of 
the common apparent horizon. The Goddard waveform 
refers to an equal-mass nonspinning binary. The simula- 
tion lasts about ~ 1516M, and the waveform has ~ 16 
cycles before merger. 

Since we will present results from these two waveforms 
it is useful to first compare them by computing the FFq. 
Although the binary parameters considered by Pretorius 
and Goddard are slightly different, we expect the wave- 
forms, especially around the merger stage, to be fairly 
close. Comparisons between (shorter) waveforms com- 
puted with moving punctures and generalized-harmonic 
gauge were reported in Ref. [1^, where the authors dis- 
cussed the different initial conditions, wave extraction 
techniques, and compared the phase, amplitude and fre- 
quency evolutions. Since the two simulations use dif- 
ferent initial conditions and last for different amounts of 
time we cut the longer Goddard waveform at roughly the 
frequency where the Pretorius waveform starts. In this 
way we compare waveforms that have the same length 
between the initial time and the time at which the wave 
amplitude reaches its maximum. In Fig. [21 we show the 
FFq as function of the total binary mass. Despite dif- 
ferences in the two simulations the FFqS are rather high. 
The waveforms differ more significantly at lower frequen- 
cies. Indeed, as the total mass decreases the FFqS also 
decrease as these early parts of the waveform contribute 
more to the signal power given LIGO's PSD. 



Any waveform extracted from a numerical simulation 
will inherit truncation errors, affecting both the wave- 
form's amplitude and phase [U, HI]- To check 
whether those differences would change the results of the 
comparisons between NR and PN waveforms, we plot 
in Fig. [5] the FFqS versus total binary mass between 
two Goddard waveforms generated from a high and a 
medium resolution run [ssj . The FFgS are extremely high 
(> 0.995). 

Based on the comparisons between high and medium 
resolution waveforms, we can estimate the FFs between 
high resolution and exact waveforms. If we have several 
simulations with different resolutions, specified by the 
mesh-spacings Xi, and Xi are sufficiently small, we can 
assume that the waveforms hi are given by 

h, ^ho + x'lhd , (4) 

where n is the convergence factor of the waveform, ho is 
the exact waveform generated from the infinite resolution 
run (xq — > 0), and hd is the leading order truncation error 
contribution to the waveform and is independent of the 
mesh spacing Xi . We find that the mismatch between the 
waveforms hi and hj, 1 — FFy , scales as 

1 - FF,, oc {x- - x-f . (5) 

In the Goddard simulations, the high and medium reso- 
lution runs have mesh-spacing ratio Xh/xm = 5/6, and 
the waveform convergence rate is n = 4 The FF 

between the high resolution and exact waveforms h^ and 
ho is given by 

FFo,, = 1 - 0.87(1 - FF,„0 , (6) 

where FFhm is the FF between the high and medium 
resolution waveforms hh and hm- That is to say, the 
mismatch between hh and ho is slightly smaller than that 
between h^ and hm, where the latter can be derived from 
the FFs shown in Fig.[2l Henceforth, we shall always use 
high- resolution waveforms. A similar calculation for Pre- 
torius' waveform gives FFq/i = 1— 0.64(1— FFfem)i though 
here Xh/xm = 2/3 and n = 2. See Fig. 6 of Ref. [2l| for a 
plot of FFfim calculated from the evolution of the Cook- 
Pfeiffer initial data ^; there FFo ranges from w 0.97 for 
M/Ms = 30 to w 0.99 for M/Ms = 100. In other words, 
the mismatch between Goddard's and Pretorius' wave- 
form shown in Fig|2] is less than the estimated mismatch 
from numerical error in the latter waveform. 

We build hybrid waveforms by stitching together the 
PN and NR waveforms computed for binary systems with 
the same parameters. At the point where we connect the 



^ The plot in Ref. |21| is for "d=16" corotating Cook-Pfciffcr initial 
data, whereas the results presented here are from "d=19" initial 
data. However, the resolutions used for both sets were the same, 
and thus the mismatches should be similar, in particular in the 
higher mass range. 
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FIG. 2: FFo between NR waveforms as a function of the 
binary total-mass M. The solid curve are generated for wave- 
forms from Pretorius and the Goddard group. The longer 
Goddard waveform is shortened such that both waveforms last 
~ 671M and contain ~ 8 cycles. The dashed curve is gen- 
erated for waveforms from the high-resolution and medium 
resolution simulations of the Goddard group. All FFs are 
evaluated using LIGO's PSD. 



two waveforms, we tune the initial time to so that the 
frequency of the PN waveform is almost the same as the 
frequency of the NR waveform (there is a subtlety trying 
to match exactly the frequencies that is discussed at the 
end of this section). The initial phase 0o is then chosen 
so that the strain of the hybrid waveform is continuous 
at the connecting point. 

In Fig.[3l we show two examples of hybrid waveforms of 
an equal-mass binary. We stitch the waveforms at points 
where effects due to the initial-data transient pulse are 
negligible. We find an amplitude difference on the order 
of ~ 10% between the Goddard waveform and the re- 
stricted PN waveform. This difference is also present in 
Pretorius' waveform, but it is somewhat compensated for 
by amplitude modulations caused by eccentricity in the 
initial data. In Ref. it was shown that PN waveforms 
with 2.5PN amplitude corrections give better agreement 
(see e.g.. Fig. 12 in Ref. [1^). However, the maximum 
amplitude errors in the waveforms are also on the or- 
der of 10% [13, ill. Since neither 2PN nor other lower 
PN order corrections to the amplitudes are closer to the 
2.5PN order, we cannot conclude that 2.5PN amplitude 
corrections best approximate the numerical waves. Thus, 
we decide to use two sets of hybrid waveforms: one con- 
structed with restricted PN waveforms, and the other 
with restricted PN waveforms rescaled by a single ampli- 
tude factor, which eliminates amplitude differences with 
the NR waveforms. We shall see that the difference be- 
tween these two cases is small for the purpose of our 
tests. 

The amplitude difference between PN and NR wave- 
forms is computed at the same connecting-point GW frc- 
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FIG. 3: We show two examples of hybrid waveforms, start- 
ing from 40Hz. The PN waveforms are generated with the 
Tpn(3.5) model, and the NR waveforms in the upper and 
lower panels are generated from Pretorius' and Goddard's 
simulations, respectively. We mark with a dot the point where 
we connect the PN and NR waveforms. 



qucncy. There is another effect which causes a jump 
in the hybrid-waveform amplitude. This is a small fre- 
quency difference between PN and NR waveforms at the 
connecting point. All our NR waveforms contain small 
eccentricities HI]. As a consequence, the frequency 
evolution Lu{t) oscillates. To reduce this effect we fol- 
low what is done in Ref. (ssj and fit the frequency to a 
monotonic quartic function. When building the hybrid 
waveform, we adjust the PN frequency to match the quar- 
tic fitted frequency (instead of the oscillatory, numerical 
frequency) at the connecting point. Since the restricted 
PN amplitude is proportional to uj^^'^{t) [see Eq. dT])], 
this slight difference between ws at the connecting point 
creates another difference between the NR and PN am- 
plitudes. Nevertheless, this difference is usually smaller 
(for Goddard's waveform) or comparable (for Pretorius') 
to the amplitude difference discussed above. 
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FIG. 4: Distribution of GW signal power. In each panel, we plot a hybrid waveform (a Tpn waveform stitched to the Goddard 
waveform) in both its original form (blue curve) and its "whitened" form (red curve) [ij. We show waveforms from six 
binary systems with total masses IOMq 20Mq , 30Mq , 40Mq , 60Mq and IOOMq . The vertical lines divide the waveforms into 
segments, where each segment contributes 10% of the total signal power. 



B. Distribution of signal power in gravitational 
waveforms 

To better understand the results of the FFs between 
hybrid waveforms, we want to compute how many sig- 
nificant GW cycles are in the LIGO frequency band. By 
significant GW cycles we mean the cycles that contribute 
most to the signal power, or to the SNR of the filtered 
signal. Since GW frequencies are scaled by the total bi- 
nary mass, the answer to this question depends on both 
the PSD and the total mass of a binary. 

In Fig. H we show the effect of the LIGO PSD on the 
distribution of signal power for several waves emitted by 
coalescing binary systems with different total masses. In 
each panel, we plot a hybrid waveform (a Tpn waveform 
stitched to the Goddard waveform) in both its original 
form and its "whitened" form [3]. The whitened wave- 
form is generated by Fourier-transforming the original 
waveform i nto the frequency domain, rescaling it by a fac- 
tor of 1/ \/ Sh{f), and then inverse- Fourier-transforming 
it back to the time domain. The reference time t = is 
the peak in the amplitude of the unwhitened waveform. 
The amplitude of a segment of the whitened waveform 
indicates the relative contribution of that segment to the 
signal power and takes into account LIGO's PSD. Both 
waveforms are plotted with arbitrary amplitudes, and the 
unwhitened one always has the larger amplitude. The 
absolute amplitude of a waveform, or equivalently the 



distance of the binary, is not relevant in these figures un- 
less the rcdshift z becomes significant. In this case the 
mass of the binary is the redshifted mass {1 + z)M. Ver- 
tical lines in each figure divide a waveform into segments, 
where each segment contributes 10% of the total signal 
power. In each plot, except for the lOMQ-binary one, we 
show all 9 vertical lines that divide the waveforms into 
10 segments. In the lOM0-binary plot we omit the early 
part of the inspiral phase that accounts for 50% of the 
signal power, as it would be too long to show. 

The absolute time-scale of a waveform increases lin- 
early with total mass M; equivalently the waveform is 
shifted toward lower frequency bands. For a M = IOM0 
binary, the long inspiral stage generates GWs with fre- 
quencies spanning the most sensitive part of the LIGO 
band, around 150Hz, while for an M = IOOM0 binary, 
only the merger signal contributes in this band. Thus, 
for low-mass binary systems, most of the contribution to 
the signal power comes from the long inspiral stage of 
the waveform, while for high-mass binary systems most 
of the contribution comes from the late inspiral, merger, 
and ring-down stages. Understanding quantitatively the 
distribution of signal power will let us deduce how many, 
and which, GW cycles are significant for the purpose of 
data analysis. We need accurate waveforms from either 
PN models or NR simulations for at least those signifi- 
cant cycles. 

From Fig. m we conclude that: 
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• For an M = IOMq binary, the last 25 inspiral cy- 
cles, plus the merger and ring-down stages of the 
waveform contribute only 50% of the signal power, 
and we need 80 cycles (not shown in the figure) of 
accurate inspiral waveforms to recover 90% of the 
signal power. For an M = 2QMq binary, the last 
23 cycles, plus the merger and ring-down stages of 
the waveform contribute > 90% of the signal power, 
and current NR simulations can produce waveforms 
of such length; 

• For an M — SOA/q binary, the last 11 inspiral cy- 
cles, plus the merger and ring-down stages of the 
waveform contribute > 90% of the signal power, 
which means that, for binary systems with to- 
tal masses higher than 30Mq, current NR simula- 
tions, e.g., the sixteen cycles obtained in Ref. [ssj . 
can provide long enough waveforms for a matched- 
filter search of binary coalescence, as also found in 
Ref. P; 

• For an M = lOOM© binary, > 90% of the signal 
power comes from the last inspiral cycle, merger 
and ring-down stages of the waveform, with two 
cycles dominating the signal power. It is thus pos- 
sible to identify this waveform as a "burst" signal. 

Similar analyses can be also done for advanced LIGO and 
VIRGO. 



C. Comparing hybrid waveforms 

We shall now compute FFqS between hybrid wave- 
forms. We fix the total mass of the equal-mass binary 
in each comparison, i.e., we do not optimize over mass 
parameters, but only on phase and time. We use the 
mismatch, defined as 1 — FFg, to measure the difference 
between waveforms and we compute them for LIGO, ad- 
vanced LIGO, and VIRGO. Note that by using FFq, we 
test the closeness among hybrid waveforms that are gen- 
erated from binary systems with the same physical pa- 
rameters; in other words, we test whether the waveforms 
are accurate enough for the purpose of parameter estima- 
tion, rather than for the sole purpose of detecting GWs. 
In the language of Ref. [l^ we arc studying the faithful- 
ness of the PN templates ^ . 

Since at late inspiral stages PN waveforms are partly 
replaced by NR waveforms, differences between hybrid 
waveforms from two PN models are smaller than those 
between pure PN waveforms. In general, the more NR 
cycles we use to generate hybrid waveforms, the less the 



^ Following Ref. Il9ll . faithful templates are templates that have 
large overlaps, say ^ 96.5%, with the expected signal maximiz- 
ing only over the initial phase and time of arrival. By contrast 
when the maximization is done also on the binary masses, the 
templates are called effectual. 



difference is expected to be between these hybrid wave- 
forms. This is evident in Figs. O [H] where we show mis- 
matches between hybrid waveforms for binary systems 
with different total masses as a function of the number 
of NR cycles n. Specifically, the mismatches are taken 
between two hybrid waveforms generated from the same 
NR waveform (from the Goddard group, taking the last 
n cycles, plus merger and ring-down) and two different 
PN waveforms generated with the same masses. 

The mismatches are lower for binary systems with 
higher total masses, since most of their signal power 
is concentrated in the late cycles close to merger (see 
Fig. [4]). Comparing results between LIGO, advanced 
LIGO and VIRGO, we see that for the same waveforms 
the mismatches are lowest when evaluated with the LIGO 
PSD, and highest when evaluated with the VIRGO PSD. 
This is due to the much broader bandwidth of VIRGO, 
especially at low frequency: the absolute sensitivity is not 
relevant; only the shape of the PSD matters. In VIRGO, 
the inspiral part of a hybrid waveform has higher weight- 
ing in its contribution to the signal power. As already 
observed at the end of Sec. [Til we can sec also that the 
difference between the Epn(3.5) and Tpn(3.5) models is 
smaller than that between the Tpn(3) and Tpn(3.5) mod- 
els. 

Figures iniin] show good agreement among hybrid wave- 
forms. In Sec. llVi as a further confirmation of what was 
found in Refs. |lol. IT]|. we shall see that PN waveforms 
from Tpn and Epn models have good agreement with the 
inspiral phase of the NR waveforms. Therefore, we argue 
that hybrid waveforms are likely to have high accuracy. 
In fact, for the late evolution of a compact binary, where 
NR waveforms are available, the PN waveforms are close 
to the NR waveforms, while for the early evolution of the 
binary, where we expect the PN approximations to work 
better, the PN waveforms (from Tpn and Epn models) 
are close to each other. Based on these observations, 
we draw the following conclusions for LIGO, advanced 
LIGO, and VIRGO data- analysis: 

• For binary systems with total mass higher than 
3OM0, the current NR simulations of equal-mass 
binary systems (16 cycles) are long enough to re- 
duce mismatches between hybrid waveforms gen- 
erated from the three PN models to below 0.5%. 
Since these FFs are achieved without optimizing 
the binary parameters, we conclude that for these 
high-mass binary systems, the small difference be- 
tween hybrid waveforms indicates low systematic 
error in parameter estimation, i.e., hybrid wave- 
forms arc faithful [l^. 

• For binary systems with total mass around 
10-20Mq, 16 cycles of NR waveforms can reduce 
the mismatch to below 3%, which is usually set as 
the maximum tolerance for data-analysis purpose 
(corresponding to ~ 10% loss in event rate). By a 
crude extrapolation of our results, we estimate that 
with 30 NR waveform cycles, the mismatch might 
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FIG. 5: We show the mismatch between hybrid waveforms as a function of the number of NR waveform cycles used to generate 
the hybrid waveforms. The LIGO PSD is used to evaluate the mismatches. In the left panel, we compare the Epn(3.5) and 
Tpn(3.5) models. In the right panel, we compare the Tpn(3) and Tpn(3.5) models. From top to bottom, the four curves 
correspond to four equal-mass binary systems, with total masses WMq, 20Mq, 30Mq, and 40Mq. The dots show mismatches 
taken between hybrid waveforms that are generated with different methods. In the left panel, we adjust the amplitude of 
restricted PN waveforms, such that they connect smoothly in amplitude to NR waveforms. In the right panel, to set the 
frequency of PN waveforms at the joining point, we use the original orbital frequency, instead of the quartic fitted one. (See 
Sec, nil for the discussion on amplitude scaling and frequency fitting). 




FIG. 6: Mismatch between hybrid waveforms as a function of the number of NR waveform cycles used to generate the hybrid 
waveforms. Following the settings of Fig. (5] we show comparisons between Epn(3.5) and Tpn(3.5), and Tpn(3) and Tpn(3.5) 
models in the left and the right panels, respectively. The solid and dashed sets of curves are generated using the PSDs of 
advanced LIGO and VIRGO. In each set, from top to bottom, the three curves correspond to three equal-mass binary systems, 
with total masses 20Mq , 30Mq , and 40Mq . 



be reduced to below 1%. 

• For binary systems with total mass lower than 
IOMq, the difference between the Tpn(3) and 
Tpn(3.5) models is substantial for Advanced LIGO 
and VIRGO. Their mismatch can be > 4% and 
> 6% respectively (not shown in the figure) . In this 
mass range, pursuing more NR waveform cycles in 
late inspiral phase docs not help much, since the 



signal power is accumulated slowly over hundreds 
of GW cycles across the detector band. Neverthe- 
less, here we give mismatches for FFqS which are 
not optimized over binary masses. For the purpose 
of detection only, optimization over binary param- 
eters leads to low enough mismatches (see also the 
end of Sec. [11]). In the language of Ref. [l^ hybrid 
waveforms for total mass lower than IOMq are ef- 
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fectual but not faithful. 



IV. MATCHING NUMERICAL WAVEFORMS 
WITH POST-NEWTONIAN TEMPLATES 

In this section, we compare the complete inspiral, 
merger and ring-down waveforms of coalescing compact 
binary systems generated from NR simulations with their 
best-match PN template waveforms. We also compare 
hybrid waveforms with PN template waveforms for lower 
total masses, focusing on the late inspiral phase pro- 
vided by the NR waveforms. We test seven families of 
PN templates that either have been used in searches 
for GWs in LIGO (see e.g., Refs. [10, SH)- or are 
promising candidates for ongoing and future searches 
with ground-based detectors. We evaluate the perfor- 
mance of PN templates by computing the FFs maxi- 
mized on phase, time and binary parameters. As we 
shall see, for the hybrid waveforms of binary systems 
with total mass M < 3QMq, both the time-domain fam- 
ilies Tpn(3.5) and Epn(3.5), which includes a superposi- 
tion of three ring-down modes, perform well, confirming 
what found in Refs. [l^, |Tl[ . The standard stationary- 
phase-approximated (SPA) template family in the fre- 
quency domain has high FFs only for binary systems with 
M < 20Mq. After investigating in detail the GW phase 
in frequency domain, and having understood why it hap- 
pens (see Sec. IIVB 2p . we introduce two modified SPA 
template families (defined in Sec. IIV B 2p for binary sys- 
tems with total mass M > SOA/q. Overall, for masses 
M > 30Mq, the Epn(3.5) template family in the time do- 
main and the two modified SPA template families in the 
frequency domain exhibit the best-match performances. 



A. Numerical waveforms and post-Newtonian 
templates 

For binary systems with total mass M > 3QMq, the 
last 8-16 cycles contribute more than 80-90% of the sig- 
nal power, thus in this case we use only the NR wave- 
forms. By contrast, for binary systems with total mass 
10 < M < 3OM0, for which the merger and ring-down 
phases of the waveforms contribute only 1-10%, we 
use the hybrid waveforms, generated by stitching Tpn 
waveforms to the Goddard NR waveforms. 

We want to emphasize that FFs computed for differ- 
ent target numerical waveforms can not directly be com- 
pared with each other. For instance, the Goddard wave- 
form is longer than the Prctorius waveform, and the FFs 
arc sometime slightly lower using the Goddard waveform. 
This is a completely artificial effect, due to the fact that 
it is much easier to tune the template parameters and 
obtain a large FF with a shorter target waveform than a 
longer one. 

We consider seven PN template families. The two 
time-domain families introduced in Sec. |TT] are: 



• Tpn(3.5) [10, [Ml: The inspiral Taylor model. 

• Epn(3.5) 0, [13, [11, Hill: The BOB model which 
includes a superposition of three quasi-normal modes 
(QNMs) of the final BH. These are labeled by three in- 
tegers {l,m,n) nil: the least damped QNM (2,2,0) and 
two overtones (2, 2, 1) and (2, 2, 2). The ring-down wave- 
form is given as: 

2 

hQNMit) = ^ A„e-(*-*-^)/"==" cos [a;22„(t - Wd) + 0„] , 

(7) 

where ujimn and Timn are the frequency and decay time of 
the QNM (I, m, n), determined by the mass Mf and spin 
0/ of the final BH. The quantities An and 4>n in Eq. ([7]) 
are the amplitude and phase of the QNM (2, 2, n). They 
are obtained by imposing the continuity of and h^, 
and their first and second time derivatives, at the time of 
matching tmatch- Besides the mass parameters, our Epn 
model contains three other physical parameters: ej, cm 
and ej. The parameter et takes into account possible dif- 
ferences between the time tend at which the EOB models 
end and the time tmatch at which the matching to ring- 
down is done. More explicitly, we set imatch = (l+et)*ond, 
and if et > 0, we extrapolate the EOB evolution, and set 
an upper limit for the et search where the extrapolation 
fails. The parameters cm and ej describe possible differ- 
ences between the values of the mass Mend = -E-cnd and 
angular momentum Send = Jcnd/M^^^ at the end of the 
EOB inspiral and the final BH mass and angular mo- 
mentum. (The end of the EOB inspiral occurs around 
the EOB light-ring.) The differences are due to the fact 
that the system has yet to release energy and angular 
momentum during the merger and ring-down phase be- 
fore settling down to the stationary BH solution. If the 
total binary mass and angular momentum at the end of 
the EOB inspiral are Mond and Jond, we set the total 
mass and angular momentum of the final stationary BH 
to be Mf = (1 - eA/)A/cnd and J/ = (1 - ej) Jond, and 
use a/ = Jf/Mf to compute uiimn and r/m„. We consider 
the current Epn model with three parameters et, cm and 
ej, as a first attempt to build a physical EOB model for 
matching coherently the inspiral, merger and ring-down 
phases. Since the e-parametcrs are related to physical 
quantities, e.g., the loss of energy during ring-down, they 
are functions of the initial physical parameters of the bi- 
nary, such as masses, spins, etc. In the near future we 
expect to be able to fix the e- values by comparing NR and 
(improved) EOB waveforms for a large range of binary 
parameters. 

We also consider five frequency-domain models, in 
which two (modified SPA models) are introduced later 
m Sec. IIVB 2[ and three are introduced here: 

• SPAc(3.5) [H]: SPAc PN model with an appropriate 
cutoff frequency /cut [10, [13 ; 
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TABLE I: FFs between hybrid waveforms [Tpn(3.5) waveform stitched to the Goddard waveform] and PN templates. In the first 
row, the two numbers in parentheses are the percentages of the signal-power contribution from the 16 inspiraling NR cycles and 
the NR merger/ring-down cycles. (The separation between inspiral and merger/ring-down is obtained using the BOB approach 
as a guide, i.e., we match the Epn(3.5) model and use the BOB light-ring position as the beginning of the merger phase.) In 
the PN-template rows, the first number in each block is the FF, and the numbers in parentheses are template parameters that 
achieve this FF. The last number in each block of the Tpn(3.5) and Bpn(3.5) models is the ending orbital frequency of the 
best- match template. For the Bpn model, the ending frequency is computed at the point of matching with the ring-down phase, 
around the BOB light ring. 



• BCV (SlJ: BCV model with an amplitude correction 
term (1— a/^/"^) and an appropriate cutoff frequency /cut- 

• BCVinipr [HI: Improved BCV model with an am- 
plitude correction term (1 — a/^/^) and an appropri- 
ate cutoff frequency /cut- We include this improved 
BCV model because Ref. [l^ found a deviation of the 
Fourier-transform amplitude from the Newtonian predic- 
tion 1^ durin g th e merger and ring-down phases (see 
Fig. 22 of Ref. [l0|). Here we shall assume n = —2/3 
in the /" power law to get the (1 — a/^/^) form of the 
amplitude correction. While it was found that the 
value of n is close to —2/3 for the ^ = 2,to = 2 wave- 
form, this value varies slightly if other multiple moments 
are included and if binary systems with different mass ra- 
tios are considered. Finally, the a parameter is expected 
to be negative, but in our actual search it can take both 
positive and negative values. 
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af/Mf 

FIG. 7: Frequencies and decay times of the least damped 
QNM 220, and two overtones 221 and 222. The scales of the 
frequency and the decay time are listed on the left and right 
sides of the plot, respectively. 



B. Discussion of fitting- factor results 

In Table U we list the FFs for hybrid target waveforms 
and three PN template families: Tpn(3.5), Epn(3.5), 
and SPAc(3.5), together with the template parameters at 
which the best match is obtained. As shown in the first 
row, in this relatively low-mass range, i.e. lOil/© < M < 
3OM0, the merger/ring-down phases of the waveforms 
contribute only a small fraction of the total signal power, 
while the last 16 inspiraling cycles of the NR waveform 
contribute a significant fraction. Therefore, confirming 
recent claims by Refs. [l3, [lli, we can conclude that 
the PN template families Tpn(3.5) and Epn(3.5) have 
good agreement with the inspiraling NR waveforms. The 
Tpn(3.5) model gives a low FF for M — 3OM0 because 



for these higher masses the merger /ring- down phases, 
which the Tpn model does not include, start contribut- 
ing to the signal power. Note that both time-domain 
templates give fairly good estimates of the mass parame- 
ters. The SPAc(3.5) template family gives FFs that drop 
substantially when the total binary mass increases from 
10A'/q to 3OM0, indicating that this template family can 
only match the early, less relativistic inspiral phase of 
the hybrid waveforms. Nevertheless, it turns out that by 
slightly modifying the SPA waveform we can match the 
NR waveforms with high FFs (see Sec. IIVB"2]) . 

In Table HIl we list the FFs for full NR waveforms 
and five PN template families: Epn(3.5), SPA^'''(3.5), 
SPA^(4), BCV, and BCVimpr, together with the tem- 
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TABLE H: FFs between NR waveforms and PN templates which include merger and ring-down phases. The upper table uses 
Pretorius' waveform, and the lower table uses Goddard's high-resolution long waveform. The first number in each block is the 
FF, and numbers in parentheses are template parameters that achieve this FF. 



plate parameters at which the best match is obtained. 
The SPA°'^*(3.5) and SPA^(4) famihes are mod ified ver- 
sions of the SPA family, defined in Sec. IIVB 21 

We shall investigate these results in more detail in the 
following sections. 

1. Effective-one-body template performances 

The Epn model is the only available time-domain 
model that explicitly includes ring-down waveforms. It 
achieves high FFs > 0.96 for all target waveforms, con- 
firming the necessity of including ring-down modes and 
proving that the inclusion of three QNMs with three 
tuning parameters et, eM and ej is sufficient for detec- 
tion. As we see in Table the values of the tuning 
parameters em and ej, where the FFs are achieved, arc 
different from their physical values. For reference, the 
Goddard numerical simulation predicts Mf ~ 0.95M 



and aj = Jf/Mj ~ 0.7 [H, and Epn(3.5) predicts 
Mond = 0.967 and aond = Jond/M^^^ = 0.796, so the two 
tuning parameters should be £m — 1.75% and ej ~ 11%. 
In our search, e.g., for M = 30Mq, ej tends to be tuned 
to its lowest possible value and et tends to take its high- 
est possible value, indicating that pushing the end of the 
Epn(3.5) inspiral to a later time gives higher FFs. 

Since the parameters em and ej depend on the QNM 
frequency and decay time, we show in Fig. [7] how uimn 
and Timn vary as functions of a/ [4^ for the three modes 
used in the Epn(3.5) model. The frequencies uJimn of the 
three modes are not really different, and grow monotoni- 
cally with increasing a/. The decay times r/„i„, although 
different for the three modes, also grow monotonically 
with increasing a/. Thus, the huge loss of angular mo- 
mentum ej, or equivalently the small final BH spin re- 
quired in the Epn(3.5) model to achieve high FFs, indi- 
cates that low ring-down frequencies and/or short decay 
times are needed for this model to match the numerical 
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FIG. 8: Frequency evolution of waveforms from the Epn(3.5) model, and the NR simulations of the Goddard group. In the left 
and right panels, we show frequency evolutions for two equal-mass binary systems with total mass 30Mq and IOOM0 . In each 
panel, there are two nearly monotonic curves and two oscillatory curves, where the former are frequency evolutions and the 
latter are binary coalescence waveforms. The solid curves (blue) are from the NR simulations, while the dashed curves (red) 
are from the Epn(3.5) model. The vertical line in each plot shows the position where the three-QNM ring-down waveform is 
attached to the FOB waveform. 



merger and ring-down waveforms. 

In Fig. [HI we show Goddard NR and Epn(3.5) wave- 
forms, as well as their frequency evolutions, for two 
equal-mass binary systems with total masses SOMq and 
IOOM0. In the low-mass case, i.e., M = SOM©, since the 
inspiral part contributes most of the SNR, the Epn(3.5) 
model fits the frequency and phase evolution of the NR 
inspiral well, with the drawback that at the joining point 
the EOB frequency is substantially higher than that of 
the NR waveform. Then, in order to fit the early ring- 
down waveform which has higher amplitude, the tuning 
parameters have to take values in Table [ll] such that the 
ring-down frequency is small enough to get close to the 
NR frequency during early ring-down stage, as indicated 
in Fig. [5] The late ring-down waveform does not con- 
tribute much to the SNR, and thus it is not too sur- 
prising that waveform optimizing the FF does not ade- 
quately represent this part of the NR waveform. In the 
higher mass case, M = IOOA/q, the Epn(3.5) model gives 
a much better, though not perfect, match to the merger 
and ring-down phases of the NR waveform, at the ex- 
pense of misrepresenting the early inspiral part. Again, 
this is not unexpected considering that in this mass range 
the merger and ring-down waveforms dominate the con- 
tribution to the SNR. 

Comparing the two cases discussed above, we can see 
that with the current procedure of matching the inspiral 
and ring-down waveforms in the EOB approach it is not 
possible to obtain a perfect match with the entire NR 
waveform. However, due to the limited detector sensi- 
tivity bandwidth, the EEs are high enough for detection. 



The large systematic error in estimating the physical pa- 
rameters will be overcome by improving the EOB match- 
ing procedure during the inspiral part, and also by fixing 
the e-parameters to physical values obtained by compar- 
ison with numerical simulations. 

Finally, in Figs. [HI UHl we show the frequency-domain 
amplitude and phase of the NR and EOB waveforms. 
Quite interestingly, we notice that the inclusion of three 
ring-down modes reproduce rather well the bump in the 
NR frequency-domain amplitude. The EOB frequency- 
domain phase also matches the NR one very well. 



2. Stationary-phase-approximated template performances 

Figures [H [10] also show the frequency-domain phases 
and amplitudes for the best-match SPAc(3.5) waveforms. 
We see that at high frequency the NR and SPAc(3.5) 
phases rise with different slopes ^. Based on this obser- 
vation we introduce two modified SPA models: 

• SPAf *(3.5): SPAc PN model with unphysical values 
of rj and an appropriate cutoff frequency /cut- The range 
of the symmetric mass-ratio 77 = toi7ti2 / (mi -t-m2)^ is ex- 
tended from its physical range ^ 0.25 to the unphysical 



^ By looking in detail at the PN terms in the SPAc (3.5) phase, we 
find that the difference in slope is largely due to the logarithmic 
term at 2.5PN order. 
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M=30Mp, M=100Mf 




Frequency (Hz) Frequency (Hz) 



FIG. 9: For M = SOMq equal-mass binary systems, we com- 
pare the phase and ampUtude of the frequency-domain wave- 
forms from the SPAc models and NR simulation (Goddard 
group). We also show the amplitude of the waveform from 
the Epn(3.5) model. 



FIG. 10: For M = 100A/q equal-mass binary systems, we 
compare the phase and amplitude of the frequency-domain 
waveforms from the SPAc models and NR simulation (God- 
dard group). We also show the amplitude of the waveform 
from the Epn(3.5) model. 



range 1. 

• SPA^(4): SPAc FN model with an ad hoc 4PN order 
term in the phase, and an appropriate cutoff frequency 
/cut- The phase of the SPA model is known up to the 
3.5PN order (see, e.g., Eq. (3.3) of Ref. [33i): 

m = 2nfto - 00 - J + ^ E a.v'^ , (8) 

where v = {nM/y/^. The PN coefficients a^s, k ~ 
0, . . . , iV, (with iV = 7 at 3.5PN order) are given by Eqs. 
(3.4a), (3.4h) of Ref. We add the following term at 
4PN order: 

as = y\ogv, (9) 

where 3^ is a parameter which we fix by imposing high 
matching performances with NR waveforms. Note that a 
constant term in as only adds a 4PN order term that is 
linear in /, which can be absorbed into the 27r/to term. 
Thus, to obtain a nontrivial effect, we need to introduce 
a logarithmic term. The coefficient y could in princi- 
ple depend on rj. We determine y by optimizing the 
FFs of equal and unequal masses. We find that in the 
equal-mass case y does not depend significantly on the 
binary total mass and is given by = 3923. The latter 
is also close to the best match value obtained for unequal 
masses. More specifically, it is within 4.5% for binary sys- 
tems of mass ratio 1712/ mi = 2. To further explore the de- 
pendence of y on rj, we need a larger sample of waveforms 



for unequal-mass binary systems As seen in Table HIl 
the two modified SPAc template families have FF> 0.97 
(except for one 0.9693) for all target waveforms, even 
though no explicit merger or ring-down phases are in- 
cluded in the waveform. The SPA^ (4) model provides 
also a really good estimation of parameters. 

In Fig. [fl] we plot Goddard NR and SPA^(4) wave- 
forms for two equal-mass binary systems with total 
masses M = 30Mq and M = lOOM©. We can clearly see 
tail-like ring-down waveforms at the end of the SPA^ (4) 
waveforms, which result from the inverse Fourier trans- 
form of frequency domain waveforms that have been cut 
at / = /cut- This well-known feature is called the Gibbs 
phenomenon. At first glance, it may appear surpris- 
ing that the often inconvenient Gibbs phenomenon 
can provide reasonable ring-down waveforms in the time 
domain. However, by looking at the spectra of these 
waveforms in the frequency domain (see the amplitudes 
in Figs. [9] & [lOl), we see that the SPA^(4) cuts off at 
the frequency /cut (obtained from the optimized FF) 
where the NR spectra also start to drop. Thus, even 
though the frequency-domain SPAc waveforms are dis- 
continuous, while the frequency-domain NR waveforms 
are continuous (being combinations of Lorentzians), the 
SPAc time-domain waveforms contain tails with frequen- 
cies and decay rates similar to the NR ring-down modes. 
We expect that the values of the cutoff frequency /cut 
at which the FFs are maximized are well-determined by 



^ Note that the auxiliary phase introduced in Eq. (239) of Ref. [43II 
also gives rise to a term in the SPA phase of the kind / logi), 
except an order of magnitude smaller than y. 
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the highest frequency of the NR waveforms, i.e. by the 
frequency of the fundamental QNM. In the next section, 
we shall show quantitative results to confirm this guess. 



3. Buonanno-Chen-Vallisneri template performances 

In Table Hi] we see that the BCV and BCVimpr fam- 
ilies give almost the same FFs for relatively low-mass 
binary systems (M = 30, 4OM0), while the BCVimpr 
family gives slightly better FFs for higher mass binary 
systems {M = 60, IOOMq). For higher-mass binary sys- 
tems, we find that the a parameter takes negative values 
with reasonable magnitude. This is because the ampli- 
tude of the NR waveforms in the frequency domain de- 
viates from the /"^^^ power law only near the merger, 
which lasts for about one GW cycle. This merger cycle 
is important only when the total mass of the binary is 
high enough (see Fig.[4|). [See also Ref. [i^ where similar 
tests have been done.] 

The BCV and BCVimpr template families give FFs 
nearly as high as those given by the SPA^ (4) family, 
but the latter has the advantage of being parametrized 
directly in terms of the physical binary parameters, and 
it gives fairly small systematic errors. 



C. Frequency-domain templates for inspiral, 
merger and ring-down 

In this section, we extend our comparisons between 
the SPAc families and NR waveforms to higher total- 
mass binary systems (40AfQ to 120AfQ) and to unequal- 
mass binary systems with mass-ratios TO2/TO1 = 1.5 and 
2. The numerical simulations for unequal-mass binary 
systems are from the Goddard group. They last for ~ 
373M and ~ 430Af , respectively, and the NR waveforms 
have ~ 4 cycles before the merger. 

In Figs, m and [13] we show the FFs for SPAf *(3.5) 
and SPA^ (4) templates, and the values of /cut that 
achieved these FFs ^. For all mass combinations (ex- 
cept for M = 4OM0 for artificial reasons) the FFs of 
SPAc(3.5) templates are higher than 0.96, and the FFs 
of SPA^ (4) templates arc higher than 0.97, confirming 
that both families of templates can be used to search for 
GWs from coalescing binary systems with equal-masses 
as large as I2OM0 and mass ratios TO2/mi = 2 and 1.5. 
Figure [13] shows that all the /cut values from our searches 



^ Note that because of the short NR waveforms for unequal-mass 
binary systems, we need to search over the starting frequency of 
templates with a coarse grid, and this causes some oscillations in 
our results. The oscillations are artificial and will be smoothed 
out in real searches. For instance, the drop of FFs at 40Mq 
for unequal-mass binary systems happens because the NR wave- 
forms are too short and begin right in the most sensitive band 
of LIGO. 



are within 10% larger than the frequency of the funda- 
mental QNM a;220 of an equal-mass binary. We have 
checked that if we fix /cut = l-07ci;22o/27r, the FFs drop 
by less than 1%. 

In Fig. 1151 we show the same information as in Fig. [7] 
except that here we draw iui„m and Tinm as functions of 
the mass-ratio 77 of a nonspinning binary. We compute 
the spin of the final BH in units of the mass of the fi- 
nal BH using the quadratic fit given by Eq. (3.17a) of 
Ref. [13: 

~ 3.35277 - 2.461r;2 . (10) 

Mf 

As Fig. ll5l shows. ^220 does not change much, confirming 
the insensitivity of the /cut on rj. 

However, in real searches we might request that the 
template family have some deviations from the wave- 
forms predicted by NR. For example, a conservative tem- 
plate bank might cover a region of /cut ranging from the 
Scliwarzschild innermost stable circular orbit (ISCO) fre- 
quency, or the innermost circular orbit (ICO) frequency 
determined by the 3PN conservative dynamics, up to a 
value slightly higher than the frequency of the fundamen- 
tal QNM. The number of templates required to cover 
the /cut dimension depends on the binary masses. We 
find that to cover the /cut dimension from the 3PN ICO 
frequency to the fundamental QNM frequency with an 
SPAJ;''*(3.5) template bank, imposing a mismatch < 0.03 
between neighboring templates, we need only two (~20) 
templates if M = 3OM0 (M = IOOM0) and 77 = 0.25. 
In the latter case, the match between templates is more 
sensitive to /cut since most signal power comes from 
the last two cycles, sweeping through a large frequency 
range, right in LIGO's most sensitive band. The num- 
ber of templates directly affects the computational power 
needed, and the false-alarm rate. Further investigations 
are needed in order to determine the most efficient way 
to search over the /cut dimension. 

For the purpose of parameter estimation, Fig. ll4l shows 
that the SPA^ (4) templates are rather faithful, giving 
reasonable estimates of the chirp mass: systematic errors 
less than about 8% in absolute value for binary systems 
with M = 4OM0 up to M = I2OM0. A difference of 
~ 8% may seem large, but the SPA^(4) templates are 
not exactly physical, and more importantly, for large- 
mass binary systems, most of the information on the 
chirp mass comes only from the last cycle of inspiral. We 
notice that when the total binary mass is higher than 
I2OM0, the FFs are relatively high (from 0.93 to 0.97), 
and the estimates of the chirp mass are still good (within 
10%). However, for binary systems with such high to- 
tal masses, the ring-down waveform dominates the SNR, 
and the SPA^ (4) template family becomes purely phe- 
nomenological. A direct ring-down search might be more 
efficient. 

All results for unequal-mass binary systems are ob- 
tained using the C22 component of ^4 (To| . which is 
the leading order quadrupole term contributing to the 
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FIG. 11: Binary coalescence waveforms from the SPA^(4) model, and the NR simulations of the Goddard group. In the left 
and right panels we show waveforms for two equal-mass binary systems with total mass 30Mq and IOOMq. The solid lines 
show the waveforms from the NR simulation, and the dashed lines give the best-matching waveforms from the SPA^(4) model. 
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FIG. 12: FFs as functions of the total binary mass. The FFs 
are computed between either the SPA^''*(3.5) or the SPA^(4) 
templates and the NR waveforms for equal-mass and unequal- 
mass binary systems. 



FIG. 13: Cutoff frequencies as functions of the total binary 
mass. We show the best-match /cut for SPA^''*(3.5) and 
SPA^(4) templates of Fig. [H The sohd black curve is the 
fundamental QNM frequency tJ22o/27r. The frequencies are in 
units of Hz. 



GW radiation. For unequal-mass binary systems, higher- 

order multipoles can also be important, and we need to ^t doing so. The Lorentzian £ is obtained as a Fourier 

test the performance of the template family directly using transform of a damped smusoid, e.g., for the fundamental 

^-4. For *4 extracted in the direction perpendicular to ^^^^ 

the binary orbit, we verified that higher-order multipoles ^ 

do not appreciably change the FFs. / gi2ir/t f ^±iuj22ot~\t\/T22o\ _ 

A natural way of improving the SPA^, models would be J ~oo ^ 

to replace the discontinuous frequency cut with a linear 2/t22o _ 2C^ if) ill] 

combination of Lorcntzians. We show here a first attempt I/tIjq + (27r/ ± a;22o)^ ~ 220U; ; 
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FIG. 14: Systematic errors of the chirp mass as functions 
of the total binary mass when SPA^ (4) templates are used. 
We show errors of the chirp masses that optimize the FFs of 
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FIG. 15; Frequencies and decay times of the least damped 
QNM 220, and two overtones 221 and 222. The scales of the 
frequency and the decay time are listed on the left and right 
sides of the plot, respectively. 



and the (inverse) Fourier transform of Eq. ([7]) reads 



E 



hQ^uif) = A„ [£+„(/) e^*" + £^,M) 



(12) 

Restricting to positive frequencies we only keep the 
^22nif) terms. In the frequency domain we attach the 
fundamental mode continuously to the SPA^ (4) wave- 
form at the ring-down frequency W220 by tuning the am- 
plitude and phase ^0 and (pQ. We denote this model 



SPA£i (note that we also need to introduce the mass- 
parameter of the final BH as a scale for ^220 and T22o)- 
Similarly, we define the SPA^a model where all three 
QNMs are combined. With the three amplitudes and 
phases as parameters, this model is similar to the spin- 
BCV template family [30| and we can optimize automat- 
ically over the 6 parameters. As an example, we compute 
the FFs between the SPA£i(4) or SPA£3(4) and the NR 
waveform of an equal-mass M = lOOM© binary. Using 
the LIGO PSD, we obtain 0.9703 and 0.9817, respec- 
tively. Those FFs are comparable to the FFs obtained 
with the simpler SPAc model, shown in Fig. [T^l It is 
known that adding more parameters increases the FFs 
but also increases the false-alarm probability. By fur- 
ther investigation and comparison with NR waveforms 
our goal is to express the phase and amplitude parame- 
ters of the Lorentzian in terms of the physical binary pa- 
rameters, relating them to the amplitudes and phases of 
the QNMs and the physics of the merger. Those parame- 
ters are somewhat similar to the e-parameters introduced 
above for the EOB model when modeling the merger and 
ring-down phases. 

We wish to emphasize that the results we presented in 
this section are preliminary, in the sense that we consid- 
ered only a few mass combinations and the NR waveforms 
of unequal-mass binary systems are quite short. Never- 
theless, these results are interesting enough to propose a 
systematic study of the efficiency of these template fam- 
ilies through Monte Carlo simulations in real data. 



V. CONCLUSIONS 

In this paper we compared NR and analytic waveforms 
emitted by nonspinning binary systems, trying to under- 
stand the performance of PN template families developed 
during the last ten years and currently used for the search 
for GWs with ground-based detectors, suggesting possi- 
ble improvements. 

We first computed FFqS (maximized only on time 
and phase) between PN template families which best 
match NR waveforms [13,1111, i-e., Tpn(3), Tpn(3.5) and 
Epn(3.5). We showed how the drop in FFqS is not simply 
determined by the accumulated phase difference between 
waveforms, but also depends on the detector's PSD and 
the binary mass. Thus, waveforms which differ even by 
one GW cycle can have FFq ~ 0.97, depending on the 
binary masses (see Fig. [T]). 

We then showed that the NR waveforms from the 
high-resolution and medium-resolution simulations of the 
Goddard group are close to each other (FFq around 0.99, 
see Fig. [2]). We also estimated that the FFq between 
high-resolution and exact NR waveforms is even higher, 
based on the numerical convergence rates of the Goddard 
simulations. 

Second, by stitching PN waveforms to NR waveforms 
we built hybrid waveforms, and computed FFqS (max- 
imized only on time and phase) between hybrid wave- 
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forms constructed with different PN models, notably 
Tpn(3), Tpn(3.5) and Epn(3.5) models. We found that 
for LIGO's detectors and equal-mass binary systems with 
total mass M > BOM©, the last 11 GW cycles plus 
merger and ring-down phases contribute > 90% of the 
signal power. This information can be used to set the 
length of NR simulations. 

The FFqS between hybrid waveforms are summarized 
in Figs. [5l [6] We found that for LIGO's detectors and 
binary systems with total mass higher than lOM©, the 
current NR simulations for equal-mass binary systems 
are long enough to reduce the differences between hybrid 
waveforms built with the PN models Tpn(3), Tpn(3.5) 
and Epn(3.5) to the level of < 3% mismatch. For GW 
detectors with broader bandwidth like advanced LIGO 
and VIRGO, longer NR simulations will be needed if the 
total binary masses M < IOMq. With the current avail- 
able length of numerical simulations, it is hard to esti- 
mate from the FFs between hybrid waveforms how long 
the simulations should be. Nevertheless, from our study 
of the distribution of signal power, we estimate that for 
M < IOA/q binary systems, at least 80 NR inspiraling 
cycles before merger are needed. 

Finally, wc evaluated FFs (maximized on binary 
masses, initial time and phase) between full NR (or hy- 
brid waveforms, depending on the total binary mass) 
and several time and frequency domain PN template 
families. For time-domain PN templates and binary 
masses lOM© < M < 20Mq, for which the merger/ring- 
down phases do not contribute significantly to the to- 
tal detector signal power, we confirm results obtained in 
Refs. [illiii, notably that Tpn(3.5) and Epn(3.5) mod- 
els have high FFs with good parameter estimation, i.e., 
they are faithful. We found that the frequency-domain 
SPA family has high FFs only for binary systems with 
M < 20Mq, for which most of the signal power comes 
from the early stages of inspiral. Furthermore, we found 
that it is possible to improve the SPA family by either 
extending it to unphysical regions of the parameter space 
(as done with BCV templates) or by introducing an ad 
hoc 4PN-order constant coefficient in the phase. Both 
modified SPA families achieve high FFs for high-mass 
binary systems with total masses 30Mq < M < 120Mq. 

For time-domain PN templates and binary masses 
M > 3OM0 , we found that if a superposition of ring- 
down modes is attached to the inspiral waveform, as nat- 
urally done in the EOB model, the FFs can increase from 
~ 0.8 to > 0.9. We tested the current Epn(3.5) template 
family obtained by attaching to the inspiral waveform 
three QNMs [3 

around the EOB light-ring. In order 
to properly take into account the energy and angular- 
momentum released during the mcrger/ring-down phases 
we introduced [Tcl two physical parameters, ej\/ and ej, 
whose dependence on the binary masses and spins will 
be determined by future comparisons between EOB and 
NR waveforms computed for different mass ratios and 
spins. We found high FFs >0.96. Due to small differ- 
ences between EOB and NR waveforms during the fi- 
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FIG. 16: The sky-average SNR for LIGO and Enhanced or 
mid LIGO detector versus total mass for an equal-mass binary 
at lOOMpc. 



nal cycles of the evolution, the best-matches are reached 
at the cost of large systematic error in the merger-ring- 
down binary parameters. Thus, the Epn(3.5) template 
family can be used for detection, but for parameter es- 
timation it needs to be improved when matching to the 
ring-down, and also during the inspiral phase. The re- 
finements can be achieved (i) by introducing deviations 
from circular motion, (ii) adding higher-order PN terms 
in the EOB dynamics, (iii) using in the EOB radiation- 
reaction equations a GW energy flux closer the the NR 
flux, (iv) designing a better match to ring-down modes, 
etc.. The goal would be to achieve dephasing between 
EOB and NR waveforms of less than a few perc ent in 
the comparable-mass case, as obtained in Ref. [48| in the 
extreme mass-ratio limit. Indeed, with more accurate nu- 
merical simulations, especially those using spectral meth- 
ods [iot . it will be possible to improve the inspiraling 
templates by introducing higher-order PN terms in the 
analytic waveforms computed by direct comparison with 
NR waveforms. 

Frequency-domain PN templates with an appropriate 
cutoff frequency /cut provide high FFs (> 0.97), even 
for large masses. This is due to oscillating tails (Gibbs 
phenomenon) produced when cutting the signal in the 
frequency domain. We tested the SPAJ;''*(3.5) and the 
SPA^(4) template families for total masses up to 120Mq, 
and three mass ratios m2/mi = 1, 1.5, and 2. We always 
get FFs > 0.96, even when using a fixed cutoff frequency, 
/cut = 1.07w22o/27r. Because of its high efficiency, faith- 
fulness, i.e., low systematic error in parameter estima- 
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tion, and simple implementation, the SPA^ (4) template 
family (or variants of it which include Lorentzians) is, to- 
gether with the EOB model, a good candidate for search- 
ing coherently for GWs from binary systems with total 
masses up to 12O-A/0. 

In Fig. [m we show the sky averaged SNRs of a sin- 
gle LIGO and Enhanced or mid LIGO [13 detector, 
for an equal-mass binary at lOOMpc. The SNR peaks 
at the total binary mass M ~ ISOM© and shows the 
importance of pushing current searches for coalescing 
binary systems to M > 100 A/q. In the mass range 
3OAf0 < M < 120Mo, the SNR drops only shghtly if 
we filter the GW signal with SPAf *(3.5) or Epn(3.5) 
instead of using NR waveforms. The difference be- 
tween Epn(3.5) and SPA°'^'(3.5) is almost indistinguish- 
able. When M > 120Mq, although the SPAf *(3.5) and 
Epn(3.5) template families give fairly good SNRs, it is 
maybe not a good choice to use them as the number of 
cycles reduces to a few. The key problem in detecting 
such GWs is how to veto triggers from non-Gaussian, 
nonstationary noise, instead of matching the effectively 
short signal. This is a general problem in searches for 
short signals in ground-based detectors. 
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APPENDIX A: COMMENT ON WAVEFORMS 
OBTAINED FROM THE ENERGY-BALANCE 
EQUATION 

In adiabatic PN models, like the Tpn model used in 
this paper, waveforms are computed under the assump- 



tion that the binary evolves through an adiabatic se- 
quence of quasi-circular orbits. More specifically, one 
sets f = and computes the orbital frequency w from the 
energy-balance equation dE{uj)/dt = where E{uj) 

is the total energy of the binary system and is the 

GW energy flux. Both E{lo) and J-{Lij) are computed 
for circular orbits and expressed as a Taylor expansion 
in oj. The adiabatic evolution ends in principle at the 
innermost circular orbit (ICO) (35j . or minimum energy 
circular orbit (MECO) [s^l, where [dE/duj) = 0. 

By rewriting the energy-balance equation, u){t) can be 
integrated directly as 



" dE{u)/du 



The RHS of Eq. (jAip can be expressed as an expansion 
in powers of lo. The expanded version is widely used in 
generating adiabatic PN waveforms [l^, [13, [Hi, [1^, it 
is used to generate the so-called Tpn template family. It 
turns out that Tpn(3) and Tpn(3.5) are quite close to the 
NR inspiraling waveforms [10, [lH ■ We wonder whether 
using the energy-balance in the form of Eq. (|Aip . i.e., 
without expanding it, might give PN waveforms closer to 
or farther from NR waveforms. In principle the adiabatic 
sequence of circular orbits described by Eq. (|Aip ends at 
the ICO, so the adiabatic model should work better until 
the ICO and start deviating (with uj going to infinity) 
from the exact result beyond it. 

In Fig. [T7] we show the NR orbital frequency uj{t) to- 
gether with the PN orbital frequency obtained by solv- 
ing the unexpanded and expanded form of the energy- 
balance equation. The frequency evolution in these two 
cases is rather different, with the orbital-frequency com- 
puted from the expanded energy-balance equation agree- 
ing much better with the NR one. When many, extremely 
accurate, GW cycles from NR will be available, it will be 
worthwhile to check whether this result is still true. 
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